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The Impact-Parameter dependent Saturation Model (IP-Sat) is a simple dipole model that incor- 
O ' porates key features of the physics of gluon saturation and matches smoothly to the perturbative 

^ , QCD dipole expression at large for a given x. It was previously shown that the model gives a 

\mmJ\ • good description of HERA data suggesting evidence for gluon saturation effects at small x. The 

, model has also been applied to proton-proton and proton-rmcleus collisions and provides the basis 

for the IP-Glasma model of initial conditions in heavy ion collisions. Here we present a reanalysis of 
available data in electron-proton collisions at small Bjorken-i, including the recently released com- 
bined data from the ZEUS and HI collaborations. We first confront the model to the high precision 
combined data for the reduced cross-section and obtain its parameters. With these parameters 
fixed, we compare model results to data for the structure function F2, the longitudinal structure 
function Fl, the charm structure function Fl'^, exclusive vector meson {J/ijj^cj) and p) production 
and Deeply Virtual Compton Scattering (DVCS). Excellent agreement is obtained for the processes 
considered at small 2; in a wide range of . Our results strongly hint at universality of the IP-Sat 
dipole amplitude and the extracted impact-parameter distribution of the proton. They also provide 
a benchmark for further refinements in studies of QCD saturation at colliders. 
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I. INTRODUCTION 



> 

(N 

. Deep inelastic scattering (DIS) and exclusive diffractive processes in electron-proton (ep) collisions such 

' as exclusive vector meson production and deeply virtual compton scattering (DVCS) at small Bjorken-x are 
I . excellent probes of the high-energy limit of QCD, in particular the parton saturation regime An effective 

field theory describing this regime is the Color Glass Condensate (CGC) [Ili. A key ingredient in particle 
production in the CGC approach is the universal dipole amplitude, the imaginary part of the quark-antiquark 
' scattering amplitude on a proton or nuclear target. The evolution of the dipole amplitude with rapidity, given a 
suitable initial condition, is computed most generally by solving the Jalilian-Marian-Iancu-McLerran-Weigert- 
Leonidov-Kovner (JIMWLK) hierarchy of equations [3]. In the large Nc limit however, the coupled JIMWLK 
equations simplify to a closed form expression for the dipole amplitude known as the Balitsky-Kovchegov (BK) 
equation [1,13. 

A common problem in QCD weak coupling approaches, to which the JIMWLK equations are not immune, is 
the impact parameter dependence of distributions; different final states will have varying degrees of sensitivity to 
large impact parameters where intrinsically non-perturbative physics may be important. The impact parameter 
dependence of the dipole amplitude in the JIMWLK and BK equations is cumbersome to implement, and 
is often in practice assumed to be uniform in impact parameter. A simple dipole model that incorporates 
the physics of saturation and models the impact parameter dependence of gluon distributions is the IP-Sat 
saturation model 0-01 . This model for the dipole amplitude, whose form can be derived at the classical level 
in the CGC contains an eikonalized gluon distribution which satisfies DGLAP evolution while explicitly 
maintaining unitarity. It also matches smoothly to the high perturbative QCD limit. The impact parameter 
dependence of the amplitude is simple to implement in the IP-Sat model and is easily generalized from DIS 
off protons to DIS off nuclei 13, Hj- It is perhaps the most attractive feature of this model and allows one 



to confront a large body of HERA data on exclusive diffractive processes which cannot otherwise be described 
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simply in saturation models-'^. 

A disadvantage of this model is that it does not have the same grounding in fundamental theory as possessed 
by the rcBK and JIMWLK frameworks. In the latter, the basis for QCD evolution is the BFKL equation, 
and higher order computations can be performed systematically in the CGC. The IP-Sat model attempts to 
approach the saturation boundary via DGLAP evolution; the eikonalization of the gluon distribution represents 
higher twist contributions that are becoming important at small x. However, computations including higher 
twist evolution are difficult to perform and enormously complicated; indeed, this difficulty was what motivated 
the CGC effective field theory approach in the first place. 

With these comments and the stated caveats, we will explore here how well the IP-Sat model fares with the 
latest combined data and attempt a detailed analysis of this data. The spirit of this computation is similar 
to the extraction of parton distributions in the collinear factorization approach. However, in contrast to the 
collinear factorization approach, the number of free parameters in the dipole approach is considerably less^; in 
the IP-Sat model we have only 4 free parameters that are varied. Recently, the HI and ZEUS collaborations 
have released new combined data for inclusive DIS , [13 with unprecedentedly high precision. Because of the 
extremely small error bars of the new combined data set, it is vital to reconfront the IP-Sat dipole model with 
these data to examine the effects of the tighter constraints on model parameters. Another advantage of the 
combined datasets is that the data are for the first time given in terms of the reduced cross-section (7^, which 
is unbiased towards any theoretical assumption in extracting the structure functions F2 and F^. 

The main purpose of this paper is to reexamine the IP-Sat model in view of recent precise data from HERA 
and to obtain its free parameters from a fit. An important application of our results will be to make predictions 



for nuclei, in particular for a future Electron-Ion Collider (EIC) jl8l - l2Cll |. Another is to apply the improved 



results for the dipole amplitude to make comparisons to data in proton-proton 21 1, proton- nucleus [22| and 



nucleus-nucleus collisions [23[. We aim the paper to be self-contained. In sections II and HI, we provide a 



concise review of the formalism for calculating the total DIS cross-section, structure functions and exclusive 
diffractive processes in the color dipole approach. In section IV, we introduce the IP-Sat dipole model. In 
section V, we present our detailed numerical analysis, and the main results; these are then confronted with 
HERA data. Section VI contains a summary of the key results of this work. 



II. TOTAL DIS CROSS-SECTION AND STRUCTURE FUNCTIONS 



In the dipole picture, the scattering of the virtual photon 7* on the proton at small x has a distinct chronology 
in light front time. First, the 7* fluctuates into a quark-antiquark pair (the so-called qq dipole) of flavor /. 
It then interacts with the proton via gluon exchanges and emissions; finally the qq pair recombines to form 
the final-state virtual photon. Because the lifetime of this qq dipole at small x is much longer than its typical 
interaction time with target, the total deeply inelastic cross-section for a given x and can be written in the 



factorized form [2J, [25 1 



al}{Q^x)^2Y, f I d'bd^ f dz\^^l]^{r,z-Q^)\^N{x,r,b), (1) 

j: J J Jo 

where Af {x,r, b) is the imaginary part of forward qq dipole- proton scattering amplitude with transverse dipole 
size r and impact parameter b (see section IV). The light cone wavefunction ^^^^ for 7* fluctuations into qq is 



^ Note that in this paper we only focus on the IP-Sat model and we do not consider the b-CGC model [T3l which is an alternative 
impact-parameter dependent saturation model that has been applied to many reactions including diffractive processes 0, [T3 ] , 
proton-proton 13], proton-nucleus [T^ . and nucleus-nucleus collisions 15] . 

^ In fairness it must be noted that the collinear approaches attempt to describe nearly all available data. Here we consider only 
data for x < 10~^. 
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computable in QED with L,T denoting the longitudinal and transverse polarizations of the virtual photon: 

|4/)(,,,;Q2)|2 ^ ^mlLj2e}{e}[K,{efr)r[z' + {l-zr]+m}[Ko{efr)r}, (2a) 
^ f 

\¥/\r,z;Q')\' = ^^J2e}{AQh'{l^zr[Mejr)r}, (2b) 



with 



j^z{l-z)Q' + mj, (3) 



where z is the fraction of the light cone momentum of the virtual photon carried by the quark, m/ is the quark 
mass, aEM is the electromagnetic fine structure constant, e/ is the electric charge of a quark with flavor /, and 
Nc denotes the number of colors. For the light quarks, the gluon density is evaluated at a; = xsj (Bjorken-x), 
while for charm quarks we take x = xbj (1 + 4m^/Q^). 

The proton structure function F2 and the longitudinal structure function can be written in terms of 7*p 
cross-section as 



F2(g^x) = ^^[ap(Q^a;)+4>(Q^x)l, (4) 



Fe{Q\x) = .^^al\Q\x). (5) 



The contribution of the charm quark to the wave functions in Eqs. ^ feeds into Eqs. ([T]) and ^ directly giving 
the charm structure function Fl"^. The reduced cross-section Ur which is a directly measurable observable is 
expressed in terms of the inclusive F2 and Fj^ to be 

Or (x, g2) = F2 (x, g2) - —-^—^Fi^ {x, q2) , 

where y = Q'^ /{sx) is the inelasticity variable and ^/s denotes the center of mass energy in ep collisions. In the 
above expression, we neglected the contribution of the Z boson which is important only at very large . 



III. DEEPLY VIRTUAL COMPTON SCATTERING AND EXCLUSIVE DIFFRACTIVE 

PROCESSES 

In the dipole picture, similar to the case of the inclusive deeply inelastic scattering cross-section, the differential 
cross-section for the exclusive diffractive process ^* + p ^ E + p with a final-state vector meson E = J/^*, (/>, p 
or a real photon i? = 7 in DVCS, can be written in terms of a convolution of the dipole amplitude and the 
overlap wave- functions of photon and the exclusive final-state particle Q, 



dt 16n 



■^T,L 



+ (6) 



with the scattering amplitude"^ 



_^7>->Bp^2i Jd^rj\z jd^b {^*e^)tx c-'l''-(i-")'-l-^AA(a;,r, b) 



(7) 



where = —t with t being the squared momentum transfer. Note that the amplitude in Eq. ^ is purely 
imaginary and obtained under assumption that the dipole-proton S-niatrix is purely real. The factor (1 -f /3^) 



^ An apparent factor of 47r difference with |3 is due to our different convention for ^t,l- 
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in Eq. ^ takes into account the missing real part of amplitude in Eq. ([7]) where 13 is the ratio of the real to 
imaginary parts of the scattering amplitude, see e.g. 

0,0,1231, 



(3 = tan 



ttA 

T 
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with A = 



■Ep 



ain(l/a;) 



(8) 



For exclusive diffractive processes one should also incorporate the skewedness effect due to the fact that the 
gluons attached to the qq can carry different light-cone fractions x, x' of proton. At NLO level, in the limit that 
x' << X << 1, the skewedness effect 12 8| can be accounted for by simply multiplying the gluon distribution 
xg{x,fi^) by a factor Rg defined via [ol. |27|, 



22T+3r(7 + 5/2) , dln\xg{x,fi'^ 
— — , with 7 = -r — ; — 

r(7 + 4) ' ' ain(l/a;) 



We note that there is uncertainty with regard to how one incorporates the skewedness correction at small x, 
and the factor Rg should be regarded as a phenomenological estimate. Nevertheless, the gluon distribution 
is mainly determined from the reduced cross-section (or structure functions) alone; the choice of Rg will only 
slightly affect the parametrization of the former. 

In the case of the DVCS for the real photon production, only the transverse component of the overlap 
wavefunction contributes. Similarly to Eq. (|2a|) . after summing over the quark helicities, we have 

t N 

= ^a,^4 + (1 _ ^)2] efK^iefr)mfK,imfr)+m,}Koiefr)Koimjr)} , (10) 

where one should sum over quark flavors / = it, d, s, c. There are several different prescription for modeling the 
vector meson wavefuctions. Following Refs. [l,!^, [2^], here we assume that in analogy to the DVCS, the vector 
meson Fock space is effectively dominated by the gg^Fock component. The overlap between the photon and the 
vector meson wave functions is then expressed as [26i] 

N, 



{^*y^)T = eje ^^^^^l^ {m)Ko{tfr)Mr, z) ~ [z^ + (1 - zf] efK^{efr)drMr. z)] , (11) 

(12) 



{^*y^)L = eje 2Qz{l - z) Ko{efr) 



ml ~ \ll 

MvMr, z) + — ^- Mr, z) 

Mvz[i — z ) 



where = {\lr)dr + 9^, is the meson mass and the effective charge is defined e/ = 2/3, 1/3, or l/\/2, for 
J/V', 4>, or p mesons respectively. In the photon case, the scalar part <j) is given by modified Bessel functions, 
whereas for vector mesons various quark models provide strong indication that a hadron at rest can be modeled 
by Gaussian fluctuations in transverse separation. We use the so-called boosted Gaussian wave-functions [26] 
which were found to provide a better description of data !9| . The boosted Gaussian wavefunction has several 
advantages over other commonly used models, namely it is more self-consistent and has the proper short- 
distance limit of z{\ — z) as m/ — >■ 0, and more importantly it is fully boost-invariant. The boosted Gaussian 



wavefunction has the following simple form 26 1 



/ "^27^2 2z(l-z)r2 mln'^\ 
0.,.(., .) = Nr.Ml - z) exp (^-^^ - + ^ j . (13) 

where the parameters Nt^l and TZ are determined by imposing the normalization condition and via a fit to 
experimentally measured leptonic decay width of the vector meson for the longitudinally polarized case Q . For 
the vector meson production, the gluon density is evaluated at x = xsj (l + My/Q"^^. 

We have followed closely here Refs. jS], where several improvements were implemented in the treatment 
of exclusive diffractive processes compared to the original paper of Kowalski and Teaney These include 
incorporating the contribution of the real part of the forward scattering amplitude and the skewedness effect 
due to the parton distributions being off-forward. We note also that in earlier studies, an important phase 
factor exp (i(l — z)r ■ A) arising from the non-forward wave- functions [29| was neglected. As shown by Kowalski, 
Motyka and Watt [9|, this extra phase factor is very important and we have included it in our analysis. In Eq. d?]). 
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the phase factors exp {—ib ■ A) and exp — z)r ■ A) effectively takes into account the size of vector mesons 
in impact-parameter and dipole transverse-size space, respectively. Even with these significant improvements, 
there is considerable room for further improvement in the study of off-forward gluon distributions at small x, 
especially in the region where higher twist effects are important. 



IV. IMPACT-PARAMETER DEPENDENT SATURATION DIPOLE MODEL: IP-SAT MODEL 

The common ingredient of the cross-sections in DIS, exclusive diffractive vector meson production and DVCS 
is the universal qq dipole-target amplitude. As seen in Eqs. (I6l7p . the impact- parameter dependence of the 
dipole amplitude is crucial for describing exclusive diffractive processes. For the total cross-section, the effect 
of the impact-parameter dependence of the dipole amplitude is not especially important and the 6-dependence 
can be effectively incorporated by treating it as a step function and adjusting the overall normalization. In this 
way, one can still find a good fit for the structure functions and total DIS cross-section. However, a consequence 
of a trivial 6-dcpendcnce leads to a pronounced dip in the ^-distribution of vector meson production at large \t\. 
This is not observed in data and can therefore be ruled out . The choice of the impact-parameter profile of 
the dipole amplitude entails non-perturbative effects coming into play that are beyond the perturbative non- 
linear BK or JIMWLK equation. Both of these small x evolution equations generate a power law Coulomb tail, 
which is not confining at large distances J30 [j32ll . A simple ^-dependence for the dipole amplitude is obtained 
by combining the Glauber-Mueller form [sl, [oj, |25j of the amplitude 



7V(x,r,6) = \ l-cxpi--^as{n'^)xg{x,fj.'^)TG{b) 
with a Gaussian impact parameter profile 



(14) 



where Bq is a dimensionful scale which, in some models, corresponds to the spatial string tension in QCD. In 
Eq. (IT4t . xg [x,iiq) is the gluon density evolved up to the scale /i with LO DGLAP gluon evolution (neglecting 
its coupling to quarks). Note that in the original IP-Sat fit, the number of flavors was taken to be 3 @, 1^. 
However, since the parameter Bq will be fixed with experimental data for exclusive J/^ production and because 
we also want to describe the charm structure function, we take Nf = 4. We then take the corresponding one 
loop running-coupling value of as with Aqcd = 0.156 GeV fixed by the experimentally measured value of as 
at the Z° mass. The contribution from bottom quarks is neglected. As in the original IP-Sat model, the scale 
fj,'^ is related to the dipole transverse size by 

^^^ = cy + fil (16) 

and the initial gluon distribution at the scale /Iq is taken to be 

xg{x,nl)^AgX-^^{l-x)^-^. (17) 

Following Refs. 0|, the parameter C is set fixed^ to C = 4. Thus the parameters Ag,Xg, fiQ and Bq are the 
only free parameters of our model which will be fixed by a fit to the reduced cross-section. At large values 
of My + Q^, we are in the color transparency regime and the main contribution to Eq. ^ comes from small 
dipole sizes. Therefore the ^-distribution at small dipole sizes can be approximately determined by the Fourier 
transform of Tc{h), 



da. 



dt 



(18) 



which is fully supported by the experimental data (see Sec. V). We will extract the values of Bq from the 
t-distribution of exclusive vector meson and DVCS data. 



^ The parameter C is correlated with other parameters of the model and its value cannot be uniquely determined via a fit. 
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FIG. 1: The parameters of the IP-Sat model, Ag,Xg,jj.Q and the corresponding y^/d.o.f obtained from the fit as a 
function of lower virtuality cut Q^in in the data bin selection. 
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FIG. 2: The /d.o.f of the fit as a function of quark masses. In the left panel, we take nic = 1.27 GeV and vary the 
light quark mass while in the right panel, we take rriu = and vary the charm quark mass. 



V. NUMERICAL RESULTS AND DISCUSSION 



We include in our fit the recently released data for the reduced cross-section (7^ from the combined analysis 
of the HI and ZEUS collaborations In earlier analyses for obtaining a fit for IP-Sat model [1, HI data 
was not included in the calculation in order to avoid introducing a normalization parameter between old HI 
and ZEUS data. The new combined data from HI and ZEUS collaborations are free from this problem and 
extremely precise. Unfortunately, experimental data for structure functions from other experiments such as 
muon DIS data from E665 ^1] and NMC [sS] have rather larger error bars and will spoil the high quality data 
points of (Tr from the combined HI and ZEUS collaborations, even though including those data will improve 
the of the fit'^. Therefore we discard the experimental data for the F2 structure function from E665 and 



^ We checked that indeed these fixed target data improve about 20 -i- 30%. 
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Data Bo/GeV^ 


fn„,d,s/GeV 


mJGeV 








xVd.o.f. 


cr,. 4 





1.27 


1.51 


2.308 


0.058 


298.89/259 = 1.15 


cr,. 4 


^ 


1.4 


1.428 


2.373 


0.052 


316.61/259 = 1.22 



TABLE I: Parameters of the initial gluon distribution in the IP-Sat model determined from fits to the reduced cross- 
section Or from the combined HI and ZEUS data in neutral current unpolarized e^p scattering 1^ for data in the range 
G [0.75, 650] GeV^ and x < 0.01. Resuhs are shown for two fixed values of the charm quark masses. 




FIG. 3: Left: The gluon structure function as a function of x for various fixed values of virtuality extracted in the 
dipole saturation model (IP-Sat), CTIO (NNLO) Q and MSTW 2008 (NNLO) [11] ■ The corresponding theoretical un- 
certainties are represented with bands between solid, dashed and dotted lines for IP-Sat, CTIO and MSTW, respectively. 
Right: The gluon structure function xg [x, At^(r)) as a function of dipole transverse size r for various fixed values of x. 



NMC collaborations. We should emphasize further that our fit to extract the free parameters of the model is 
performed for the reduced inclusive DIS cross-section alone. With the extracted parameters, we then confront 
model results for F2,Fl and i^l'^ with the HERA data. 

In the calculation of x^, the statistical and systematic experimental uncertainties are added in quadrature. 
Since new combined data points have error bars as small as ~ 1%, in extracting x^, we have to evaluate all 
the theoretical observables with at least one or two orders of magnitude higher precision. This causes the 
minimization algorithms to be very slow, a problem which can be overcome by implementing a parallelization 
of the numerical code. 

As noted previously, we have altogether 4 free parameters which we fix via a fit to experimental data. 
Following Refs. [1, , we do not allow the parameter Bq to vary in addition to other parameters in the 
minimization algorithm, but we adjust it iteratively in order to obtain a good description of the t-dependence 
of exclusive diffractive processes. We only consider data within the kinematic bin with x < 10~^ and Q^i^ < 
(5^[GeV^] < 650. Although changing the upper limit cut on will not greatly influence our fit, imposing 
a lower cut on (5^ji,j significantly improves the fit. Note that in different studies, different values for a lower 
cutoff of the virtuality have been taken [1, [s, 35, 3^. Obviously, our model which is based on weak coupling 
dynamics is unreliable at very low virtualities on the order of or below the confinement scale. Moreover, the 
DGLAP evolution starts at the scale /Xq, therefore at < /ig we are sensitive to the parametrization of the 
gluon distribution at the initial scale. In Fig. [1] we show the sensitivity of our parameters Ag,Xg, and /ip 
obtained from a fit with different data bin selections corresponding to different lower cutoffs Q^in- We also 
show the corresponding x^/d-o.f as a function of Q^i^- There are two main features in Fig. [T] which guide us 
to choose the preferred value of Q^jn- Namely, /d.o.f tends to decrease by increasing and then flattens 
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FIG. 4: Left: The saturation scale extracted from the IP-Sat model with the parameter sets given in table |I] and the 
old parameterization set from Ref. Q labeled by IP-Sat (old) as a function of Xjx at various impact-parameter h. Right: 
The saturation scale as a function of the impact-parameter h for various fixed values of a::. In both panels, the lower and 
upper curves in the band correspond to the results obtained with the parameter sets given in table U with charm mass 
rua = 1.27, 1.4 GeV, respectively. 



at about ~ 0.75 GeV^. Moreover, the values of all parameters obtained from the fit become stable and 

saturate at ~ 0.75 GeV^. Therefore in the final fit shown in table U we take HERA combined data for cjr 

within £ [0.75, 650] GeV^ and x < 0.01. Altogether we include 263 experimental data points for cr^ from the 
combined HI and ZEUS data in neutral current unpolarized e^p scattering 

In Fig. [21 we show the effect of light and charm quark masses in our fit. Remarkably, the data prefer very light 
current quark masses of order a few MeV. This seems to be natural given the fact that the saturation scale here 
is large relative to Aqcd; our fit should therefore be less sensitive to non-perturbative infrared dynamics. To 
show more clearly the effect of light quark masses in our fit, we allow the light quark mass to„ (to„ — ~ mg) 
to vary as a free parameter and obtain its preferred value from the minimization. We find to„ « 10~^ GeV 
from the fit. It is seen however in Fig. [5] (left panel) that the results of the fit and minimization is stable with to„ 
of order or smaller than few MeV. In earlier analyses of the IP-Sat model, heavier quark masses of m,j = 50, 140 
MeV were required for a good fit [1, in sharp contrast to our conclusions from the new combined Hl-fZEUS 
data. 

The charm mass dependence of our fit is shown in Fig. [2] (right panel). It is seen that the minimization is 
more sensitive to light quark masses than charm mass. This is understandable given that charm data will not 
greatly influence our fit even if we include cr^'^ data in our calculation. In table HI we show our fits with two 
different charm masses rric = 1.27 and 1.4 GeV. Note that the parameters of the IP-Sat model obtained here 
from the new combined data are different from the earlier studies. If we take the old IP-Sat parameters with 
the same quark masses as in Ref. \^ and confront it with new data within £ [0.75,650] GeV^, we obtain 

/d.o.f — 825.86/259 = 3.19 which is not a good fit according to the hypothesis-testing criterion. The earlier 
IP-Sat fits should be therefore superseded by the fit given in table HI 

In Fig. [3] (left panel), we compare the gluon structure function at various fixed values of virtuality obtained 
from the dipole saturation model (IP-Sat) and the leading twist collinear factorization approach with NNLO 



DGLAP evolution, namely CTIO [371 and MSTW 2008 [38|. The bands for CTIO and MSTW correspond to 
uncertainties in obtaining a fit from global data analysis, while in the IP-Sat model the uncertainties are mainly 
due to our freedom to choose different values for the charm quark mass in the range rric = 1.27-^1.4 GeV. At large 
virtualities, we are in the color transparency region, the saturation effects become irrelevant and our approach 
approximately matches the standard perturbative formalism. The small differences seen at high virtualities are 
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10"' 10"^ 10"' 10"^ X 

FIG. 5: Results for the structure function F2{x, Q^) as a function of x for various values of . In order to separate data 
for each from the others, the data and model results represented by the lines are multiplied by a factor 2\ with i 
given on plot. We used the parameters set of the IP-Sat model given in tableUwith rUc = 1.27 GeV. The experimental 
data are from combined HI and ZEUS collaborations 1161 1. 



mainly due to the fact that we used LO DGLAP evolution without including quark degrees of freedom, while 
quark evolution contributions were included in the perturbative leading twist results shown in Fig. [3l At low 
virtualities and low x, we are in the saturation regime and we observe our gluon distributions to be significantly 
different from those obtained from the leading twist perturbative computations and significantly more stable 
than these, especially the MSTW fits. One may conclude that the higher twist contributions significantly 
influence the extraction of the gluon distribution in this regime. In Fig. [3] (right panel), we show the gluon 
structure function xg {x,iJ?{r)^ in the IP-Sat model as a function of dipole transverse-size r for various fixed 
values of x. Note that large dipole transverse-size corresponds to low virtuality via Eq. ()16|) : consequently, the 
gluon distribution is reduced. 

Following Refs. [1, @], we define the saturation scale — with being the saturation radius, as a 

scale where the dipole scattering amplitude has a value Af{x, r^, 6) = (1 — exp(— 1/2)) = 0.4. In Fig. 2] (left), 
we compare the saturation scale extracted from the IP-Sat model with old parametrization of Ref. [9J and the 
parametrization obtained in this paper as a function of 1 /x at various impact-parameter b. It is seen that the 
saturation scale extracted from old and the new combined data from HERA are consistent with each other. In 
Fig. m (right), we show the saturation scale as a function of the impact-parameter b for different fixed values of x 
obtained with parameter sets given in table HI The band in Fig. |4] shows the theoretical uncertainties associated 
with our freedom to choose different values for the charm quark mass within 1.27 ~ 1.4 GeV. It is generally 
seen that the saturation scale plotted as a function oi 1/x grows faster for more central collisions (b « 0). 
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FIG. 6: IP-Sat model results for the charm structme function F2'^{x, Q^) as a function of x for various values of 
compared to HERA data. The solid and dashed lines are obtained with two different parameter sets given in table 
U] corresponding to charm mass nic = 1.27, 1.4 GeV respectively. The experimental data points are from the recently 
released combined data sets from the HI and ZEUS collaborations !17j assuming that a^'^ ~ F2'^ (see text). 



Moreover, the saturation scale at different impact-parameter can be significantly different, as large as one order 
of magnitude. This non-trivial behavior clearly show the importance of the impact-parameter dependence of 
the saturation scale. 

With the parameters given in table U extracted from the x-squared fit to the reduced inclusive DIS cross- 
section, we now compute the structure functions F2(x,Q'^), the longitudinal structure function Fl{x,Q'^) and 
the charm structure function F|'^(a;, Q^) in the IP-Sat model (using Eqs. (1415114^ ) and compare to the combined 
HERA data sets for the same. As emphasized previously, these F2, Fl and F^"^ experimental data were not 
included in our fit and are therefore a prediction of the model based on the fit to the reduced inclusive DIS 
cross-section alone. It is seen in Figs. (|5I6I7I) that the IP-Sat model resuhs are in good agreement with structure 
functions data for a wide range of kinematics for e [0.75,650] GeV^ and x < 0.01. For the F2 and F2'' we 
confronted model results with new Hl-f ZEUS combined data while for Fl{x, Q^) we compared with the existing 
data from HERA. 

The agreement of the IP-Sat model with data is most striking for the case of structure functions F2 and Fj'^ 
where the experimental data are more precise. This agreement can also be taken as a non-trivial consistency 
check of the model. For the structure functions F2, the agreement with data extends to x > 0.01 and only begins 
to fail closer to a; = 0.1. In Fig. ([5]), we show the F2'' data from the Hl+ZEUS combined analysis assuming that 
the contribution of Fj^'^ to the reduced cross-section, originating from the exchange of longitudinally polarized 
photons, is negligible in the kinematic range considered here. This may give rise to less than a few per cent error 
to the data (not shown in Fig. (jH])). It is seen in Fig. ((T]) that there are some irregularities in the experimental 
data points for F[^{x,Q^) [39j which await more precise measurements. In Figs. (I5I6I7I) we extend our model 
results beyond the kinematics of existing data, as predictions for future DIS experiments. 

We now turn our attention to results for exclusive diffractive processes at HERA. We first focus on the t- 
distribution of exclusive vector meson production and Deeply Virtual Compton Scattering (DVCS). We fixed 
the width of proton impact-parameter profile Bq in Eq. (jl5l) via a fit to the slope of the ^-distribution of the 
J/vP mesons and we found Bq — 4GeV~^. This value is in accordance with Ref. 9] since we used the same 
data set to fix the value of Bq- In Fig. |8] we show a compilation of extracted values of the slope Bd from a fit 
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FIG. 7: IP-Sat model results for the longitudinal structure function Fl{x,Q^) as a function of s for various values of 
compared to HERA data. We used the parameters set of the IP-Sat model given in table U with = 1.27 GeV. The 
experimental data are from [39l.l4o|. 



of the form da/d\t\ cx °* for exclusive vector-mesons electroproduction and DVCS. Note that" for very large 
+ My or small dipole-size, we have Bd ~ Bq, see Eq. ([T5|) . As it is seen in Fig. [51 the experimental errors 
for the values of Bd are rather large. This leads to some uncertainties in extracting the value of the parameter 
Bg- We estimated that the uncertainties of the value of Bq is about 0.4GeV~^. 

In Fig. [51 only the line labeled by J/'J' is a fit, the other lines corresponding to light mesons and the DVCS 
are results from the model. The extracted values of Bd from theory at the same + My are larger for lighter 
vector mesons than those for J/'^ in accordance with data. This is simply because due to the convolution 
of the photon wavefunction with the vector meson wavefunction, the typical dipole sizes for light and heavy 
mesons are different and mainly controlled by 1/e/ with e/ given in Eq. ((Sj. Therefore, at a fixed virtuality, the 
typical dipole size is bigger for lighter vector meson and consequently the validity of the asymptotic expression 
in Eq. (jlSp is postponed to a higher virtuality. The recent data from ZEUS for meson production 4l| which 
extend to higher (shown also in Fig. [8]) , in accordance with the underlying assumptions of our model, indicate 
that indeed at large + My the value of Bd tends to saturate to a universal value mainly determined by the 
interaction area and shape of proton. 



In general, Bd ^ Bq because the Fourier transform of the ^-distribution in the IP-Sat exclusive VM amplitude is not a simple 
exponential in \t\. 
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FIG. 8: A compilation of the value of the slope Bd of f-distribution of exclusive vector-meson electroproduction and 
DVCS processes, as a function of + My. The compilation of experimental data is from [4l| . 
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FIG. 9: Differential vector meson cross-sections for J/^, p and (/> as a function of \t\. Data for a given W with varying 
are compared to results from the IP-Sat model using two parameter sets given in table [II wi th rUc — 1.27 GeV (solid 
line) and nic = 1.4 GeV (dashed line). The data are from the HI and ZEUS collaborations [4ll446|]. 



To further test the IP-Sat dipole model, in Figs. (I9ll0llip we confront experimental data from HI and ZEUS 
with model predictions for the \t\, and W^-dependence of the vector mesons J/^, (/> ^-nd p production in 
various kinematics. Note that none of the data sets in those figures are used into our fits and the results of the 
model are fixed by Bq alone. For the total cross-section, we performed the integral over \t\ up IGeV^. It is 
generally seen that the agreement between our results and data is excellent. In Figs. 191101 and [TT] we show our 
results calculated with two dipole parameter sets given in table U corresponding to charm mass rric = 1.27 GeV 
(solid lines) and rric — lA GeV (dashed lines). Note that both set of parameters give very similar results. The 
J/vP production cross-section is more sensitive to the charm quark mass at low compared to light vector 
meson production cross-sections and the DVCS. This is due to the fact that the scale in the integrand of the 
cross-section is set by e/ defined in Eq. ([3]), and only at low virtualities for < does the cross-section 
become sensitive to quark mass. It is seen from Fig. [TT]that the W^-dependence of the cross-section follows a 



13 



a 



10' 



y'p J/v P 



YP^PP 



YP^PP 



10" 



!» W = 90 GeV ; 




^Hi : 




+ ZEUS 







10' 



10' 

Q2+Mj/^(GeV^) 




10' 10 

Q^+Ml (GeV^) 



Q2+M2 (GeV^) 



(GeV2) 



FIG. 10: Total vector meson cross-section a for J/'if, (j) and p as a function of + My compared to results from 
the IP-Sat model using the two parameter sets given in table |T] with mc = 1.27 GeV (solid line) and rric = 1.4 GeV 
(dashed line). The ZEUS data points are scaled to HI Q'^-bins using the dependence measured by ZEUS of the form 
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FIG. 11: Total vector meson cross-section as a function of W compared to results from the IP-Sat model using the two 
parameter sets given in table U with m-c = 1.27 GeV (solid line) and rric = 1.4 GeV (dashed line). As in Fig. 1101 the 
experimental data are from HI and ZEUS collaborations [lll - G^ . 



power-law behavior of the form a oc W . Our extracted values of 5 at difFerent kinematics for various vector 



mesons are also in perfect agreement with experimental data from HI and ZEUS 4l|-|46j (for brevity, we do not 
show the comparison here) . 

In Fig. [T^ we show the ratio of the longitudinal to the transverse cross-sections R = Ohlf^T for J/^, and 
p production as a function of . The ratio i? provides vital information about the effective averaged dipole 
transverse-size which participates in the interactions. The transverse component of the cross-section is generally 
dominated by larger dipole sizes than the longitudinal component. This leads to a relative increase of the t-slope 
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FIG. 12: The ratio R = aj^joT as a function of . The solid and dashed lines are results from the IP-Sat model with 
charm mass rric = 1.27 and 1.4 GeV respectively. The experimental data are from HI and ZEUS collaborations [4ll - |46l ]. 



of the transverse component compared to longitudinal one, i.e. But > Bjjl seen also in data. Therefore one 
expects the ratio R = (Jl/<^t to increase with Q^, as can be seen in Fig. [T^l Note that the theory curves in 
Fig. [T2]were calculated at a fixed W = 90 GeV while the experimental data are for a range of M^-bin shown in 
the plot. However, the W dependence of the ratio is very weak for a fixed 
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In Fig. [m we compare the IP-Sat results for DVCS with the experimental data from HI; in the left panel, 
we show the the t-distribution at various values and fixed W = 82 GeV, in the middle panel we show cr as a 
function of with fixed W, and in right panel the total cross-section cr as a function of W with different fixed 
values of Q^. It is generally seen that the agreement with data is remarkably good. We do not compare with the 
ZEUS data, but HI and ZEUS data for DVCS are consistent with each other (43,0] and our description of ZEUS 
data is of the same quality as for HI data. The VF-dependence of the cross-section in Fig. [13] can be described 
by a power-law behavior where S « 0.76 ± 0.02 at = 8GeV^ and S « 0.83 ± 0.02 at = 25GeV^ 
extracted from the IP-Sat model is in agreement with the experimental value 6 — 0.77 ± 0.23 ± 0.19 471. l48l| . 



VI. SUMMARY 



In this paper, we confronted the IP-Sat dipole model with the new combined data from HERA which are 
significantly more precise than previous data sets. The IP-Sat model has very few free parameters and they 
are fixed from a fit to the reduced cross-section (7^ at small-z. Model results are then compared to F2, F2'^, 
and data from exclusive diffractive processes such as vector meson production and the DVCS with the available 
data from HERA including the new combined data for crj:'^. Overall, the model provides an excellent description 
of data in the range G [0.75, 650] GeV^ and x < 0.01. The quality of this description relies on a (5^j„ of 0.75 
GeV^; one observes a steady increase in x-squared with decreasing values of Qmin- Other key features of this 
novel fit are the lower values for the light quark masses, and the four flavor running of the coupling. 

In addition to the inclusive data, we showed that most features of HERA data on exclusive processes are 
reproduced by the IP-Sat model. The VF-dependence of the cross-section for vector meson production and the 
DVCS in Figs. ()ll|13p shows that the energy-dependence of the dipole amplitude is correctly modeled, while the 
|i|-dependence in Figs. (|9ll3p suggest the shape of the proton is correctly implemented in our model. Finally, 
the Q'^-dependence of the cross-section in Figs. (jl0|13|) indicates that the virtuality dependence of the gluon 
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FIG. 13: Left: DifTerential DVCS cross-section as function of \t\. Middle: Total DVCS cross-section as a function of 
Q^. The solid and dashed lines are results from the IP-Sat model with charm mass rric — 1.27 and 1.4 GeV respectively. 
Right: Total DVCS cross-section as a function of W. The experimental data are from HI and ZEUS collaborations 



distribution and its evolution with the dipole transverse-size are correctly modeled. Moreover, our model also 
reproduces the experimental data for the ratio R — (Tl/'^t shown in Fig. 1121 This is another piece of evidence in 
favor of the dipole picture and the fact that the dipole transverse size sets the relevant length scale in interactions 
for both transverse and longitudinally polarized virtual photons. We stress that the t-distribution of all three 
vector mesons {J/tp, 4>i p) ^ well as DVCS can be correctly reproduced by fixing only one parameter Bq despite 
the fact that the vector meson and DVCS wave functions are very different. This strongly hints at universality 
of the extracted impact-parameter distribution of the proton. 

Unintegrated gluon distributions can be extracted from the dipole amplitude in the IP-Sat model, and in a k± 
factorized formalism have been compared to data from RHIC and the LHC in proton-proton and proton-nucleus 



collisions (2l|, \2^. In addition, the IP-Sat model is the basis of the IP-Glasma model (23| of initial conditions 



in heavy ion collisions. However, the parameters employed in these studies were determined from data from HI 
and ZEUS predating the combined data sets for the proton. It will be interesting to see what the impact of the 
new fits are on final state observables in proton-proton, proton-nucleus and nucleus- nucleus collisions. 
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